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A nonperturbative theory of multiphonon anharmonic transitions between energy levels of a local 
mode is presented. It is shown that the rate of transitions rearranges near the critical level number 
n cr : at smaller n the process slows down, while at larger n it accelerates in time, causing a jump-like 
loss of energy followed by the generation of phonon bursts. Depending on parameters, phonons are 
emitted in pairs, triplets etc. 



PACS: 71.23.A, 71.55.-i 



The nonlinear dynamics of strong vibrational excitations in crystals has attracted remarkable attention during the 
last years. The opportunities of applying novel numerical techniques play an essential role in the investigation of 
nonlinear discrete systems, allowing for the observation of different phenomena. One of the effects, the existence 
of localized vibrations in a pure anharmonic lattice (so-called intrinsic local modes (ILM) or self-localized solitons 
[ (SLS))ErEl, is of particular interest, as it links local lattice dynamics to the physics of solitons and underlines, in 
£3 ■ general, the importance of strongly excited local modes. So far the research of strong anharmonic effects in lattice 
dynamics as well as molecular dynamics simulations of a nonlinear lattice have been mainly carried out within the 
1 . framework of classical mechanics. However, quantum effects are important here as they lead to extra mechanisms in 
the decay of vibrational excitations. This holds especially for the modes with frequencies above the phonon spectrum 
because they, being singly excited, do not decay at all in the classical limit - it is because all their harmonics are out 
' of resonance with the linear spectrum which-,makes it impossible to transfer energy to phononsEl. Quantum decay of 
the ILMs has been discussed by Ovchinnikovtl in the framework of standard perturbation theory. However, this theory 
may not be applicable here because ILMs exist only when their amplitudes (energy) are large enough; therefore, one 
cannot assume that anharmonic interactions are weak. To understand how relaxation proceeds when the energy of the 
i mode is large, one should develop a nonperturbative theory of multiphonon transitions between the quantum levels 
f*"*. ' due to anharmnnic interaction with the phonon continuum. 

Recently, irou we proposed a theory of anharmonic damping of a strongly excited mode, supposing that the-latter 
can be considered classically (phonons in this theory are described quantum- mechanically) . The idea of Ref .Lru was 
to consider the effect of the strongly excited mode upon the phonon operators. The goal was to find a transformation 
of these operators in time, and, from that, to obtain the probability of creation of a phonon. Adding all the created 
phonons allows one to find the energy loss in the mode. 

The realization of the method in the case of quadratic interaction of the mode with respect to phonon coordinate 
operators Xj is straightforward: the phonon Hamiltonian with the anharmonic interaction Viw (xi(t))xiXi' , where 
I i xi(t) = A cos (ujit) is a classical time-dependent coordinate of the local mode, can be diagonalized exactly for an 
arbitrary amplitude A of the mode. This is achieved by the transformation bj(t) = £w [/%■(£) a* + Vij(t)af] where di 
and df, are initial destruction and creation operators of phonons (for explicit expressions for fiij and z/y, sect). The 
, initial zero-point state, |0), is not the zeroth state of the operators bj. This means that phonons are generated in a 
lattice. The number of created phonons j equals nj — \ | 2 . For times t uj^ 1 the transformation simplifies^!: 
~ 5ij and the rate of energy relaxation becomes 



E,~^2nu t d\vi\ 2 /dt. (1) 



Inl it was found that for a large energy mode, the relaxation rate diminishes with energy: — Ei ~ A 2 , i.e. 
relaxation accelerates with time. This result contradicts standard perturbation theory. Yet, the interpretation of 
this surprising result is straightforward - the decay is of purely quantum mechanical origin and, therefore, it should 
disappear in the classical limit Ei — > oo. It was also shown that the acceleration of the process continues down to 
an energy E cr (which typically corresponds to Ei ~ 10 — 30 vibrational quanta, or A ~ 0.2A — 0.5A ) when the rate 
becomes high (of the order of Tiojf) and a sharp relaxation jump takes place accompanied by generation of a burst of 
phonons. The latter are emitted in pairs. After the first jump, several other jumps may occur. For small energies 
relaxation slows down exponentially in agreement with perturbation theory. 
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The theory of RefBEl considers the relaxation of the mode as a continuous process where the energy of a classical 
subsystem is diminishing. In this communication, the effects of the quantum origin of the states of the mode and of 
the transitions between its levels will be considered. Here we restrict ourselves to quadratic terms in the expansion 
of the potential energy with respect to phonon operators: U = Uq{x{) + J2 i Uu{xi)xi + Y^w U2W (xi)xiXi> /2; the 
potential of the local mode, Uq(xi), is supposed to be anharmonic, and so are the factors Uu(x{) and U2u>(x{): they 
depend on the displacement operator xi due to higher order anharmonicities. To clarify the range of the applicability 
of the approximation, we observe that all anharmonic terms become comparable to each other for shifts of an atomic 
distance of typically lAl. E.g. the harmonic V^xf and anharmonic (V^ /3 + xi/12)xf terms in the expansion of 
Uq(xi) are of the same order of magnitude for xi ~ 1A, while the terms Vfr XiEi' an d (V^f) + V l < ^)xi/2)xiXiXi' in the 
expansion of Uiu 1 (xi)xixi' are of the same order of magnitude for xi ~ 0.3A Below we take into account these terms 
while we consider just local modes with the amplitude A ~ 0.3A or less. For such amplitudes the terms ~ XiX^x^ 
are already much smaller and are neglected. 

Within this approximation, the interaction of the local mode with phonons has the form 

H^t = y]Vui{n)xj + - ^2[Viu'(n)ai +H.c.]x i x v , (2) 

i a' 

where n = af&i, af and a/ are the creation and destruction operators; terms ~ as well as ~ af , k = 2, 3, ... are 
neglected; they describe transitions between the levels n and n — k and could be considered analogously. Below the 
operators Viu(h) and Viw(n) stemming from Uu and Via 1 respectively, are applied only to the initial number state 
\n) of the local mode. Therefore the V(h)— operators can be replaced by C-functions of n. In this approximation, 
using shifted phonon operators x[ = ii + Xoi, with Xqi — Viulo^ 2 , the interaction part of the hamiltonian reads 

H' int = Vux't + 5 J2 Vm ' ^1 + a o^4. 

i ii' 

where % = Vu(af + oj) + V Ui (n -n),% = - £\, Viii'Viu^ 2 (all F-factors depend on n). 

Let us present the operators of the phonon coordinates in the form x'^t) = (h/2uji)[gi(t)a i + H.c.}. Here the 
operators gi(t) account for the time dependence of x'^t). This dependence is determined by the equations of motion 

9i + ^i9i = -Ui{Vu + Y Vui>gi>{ai + af)) (3) 

i 1 

and by the initial conditions g%{t) = e~' iu i t + i( Pi ^ 1 _ > q, fa is the random phase (the final expression should be averaged 
over fa). The operators cji do not commute with ai + af; the order of the operators in right-hand side of Eq. (3) 
depends on whether one considers the operators §i or the operators gf '; the opposite order describes gf. The integral 
form of the equations (3) is the following: 

&(t) = e-*"'-H* + / dt' sin (ui(t-t'))[V H (t') + ^(WZ'fa'Wi&itf) +&t (*'))]• 
Jo v 

Here (eiVii-i') = Vm> = J2mm' ^imVi mm 'ei'm' , Tn is the index of Cartesian coordinates of the atoms in the lattice, 

— 1/2 

e jm = &im^>i , &im are projections of the reduced displacements x m to the normal coordinates. 

We assume that relaxation of the mode is slow as compared to its frequency. Then describing the decay of the 
n-th level caused by the n — > n — 1 transition, one can neglect all other levels of the mode. This means that one 
can use the rotating wave approximation when only terms with creation and destruction operators, being arranged 
alternatingly in the iteration series for the cji, are taken into account. Let us consider in this approximation the large 
time asymptotics of the two-time correlation function (n|x^(t + r)x^(t)|n) for f > r ~ u^. Here \n) denotes the nth 
state of the mode and the ground state of phonons. Only the matrix elements (n\gi\n) and (n|<7j|n — 1) contribute 
to this asymptotics (|n — 1) = n~ x / 2 ai\n)). The terms ~ Vm(h — n) of the expression for gt(t) give zero; the terms 
~ Vu(ai + at) of this expression oscillate fast and, therefore, give small contribution which can be neglected (note, 
that the basic time-dependence of the operators ai and a/~ is given by the factors e~ tuJlt and e luJlt respectively). In 
this approximation (ra|£<(< + T)x' t {t)\n) w (h/2uj l )[\^(t)\ 2 e- tuJ ' T + \vi{t)\ 2 e^], t > r ~ uj^ 1 , where 

W (t) « e*** - (i/2) f dt'^(e i V i e i Oe i ^ +Wi ' 3t '^(i / )(n- l\ai(t')\n), 
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Vi (t) « (i/2) / < ft , y;(e i V r I e 4 Oc- <(wi+w< ' )t 'w'(*')<»|a^(* , )|n - 1) 



(/ii ~ (n|^|n), i/j ~ 1); here it was taken into account that (n\aj~ (t')\n — l)(n— l\ai(t)\n) = (n|a ; + (t')a/(i)|n)). 

This asymptotics corresponds to the following transformation of phonon operators : 6j ~ pi&i + af . To find this 
transformation explicitly we present the above expression for v% (t) in the form 

!/*(*)« (t/2) / dt'e-^'imWVinKWW-l). (4) 
Jo 

Here D{t) — J2i e — * u, **ef//i (t) satisfies the equation 

/>oo />oo 

D(t)=D (t)+ dt' dt"G(t-t')ViG*(t' -t")ViF(t' -t")D(t"), (5) 



A)0) = Ei eie-^+^S G mm ,{t) = -(i/2)6(i) £ 2 e^e^e"^* is the Green's function of phonons, 0(f) the Heavi- 
side step-function, and 

F(t) ~ <n|of a,(t)|n>. (6) 

Correlation functions of this type are well known in the theory of optical transitions (see, e.g.0): F(t) — ne~ lu)ltJr ^ t \ 
where f(t) = —Ft + 'Yl li {V^ i /2T%Lo^)(e luJit — 1); the sum accounts for phonon transitions caused by a shift in their 
equilibrium positions in connection with the n — > n — 1 transition in the mode, T = —Ei/hCui is the decay constant, 
which should be determined self-consistently (if the relaxation is slow as compared to characteristic rates of the phonon 
Green's functions, then one can neglect the term Ft in this expression); huii = hui — J2i Viu/^uf is the renormalized 
level spacing of the local mode. Equation (5) can be easily solved by means of a half-axis Fourier transform. One gets 

D{uj) = R(lu)D q (lu): R(lo) = [I - G(lu)vG*(lui - w)i)] _1 fl Q (u) (7) 

(v = Viy/n). Here the functions Do(u>), G(cu) and G*(ui) are the half-axis Fourier transforms of the functions Do(t), 
G(t) and G*(t)eX*\ 

Inserting Eqs. (4) - (7) into (1) and averaging over <pi (leading to the replacement of e*(* i_< ^') by Sw), one gets the 
following formula for the damping rate of the nth level of the mode: 

1 f°° 

T(n)~— dwSp[R(u))vp(w)R*(w)vp(u>i-u)]. (8) 
2tt Jo 

Here p(u>) = ImG{oj) and p(w) = ImG(uj) are the usual and renormalized phonon density functions, respectively. 
The damping rate of the energy is Ei >n = huji n T(n) n A comparison of this result with an analogous formula, obtained 
within the classical description of the local modejD shows their coincidence in the limit f(t) = (G(u>) = G(u))). 
Such a difference should be expected: the function f(t) accounts for phonon transitions, assisting the destruction of 
a quantum of the local mode, thus being of quantum mechanical origin. 

For moderately large n, T(n) ~ w(n) _1 , i.e. the relaxation accelerates with time, which is analogous to —E% ~ A~ 2 
obtained in the classical description of the local modal The acceleration continues down to the level n cr , when the 
resolvent in R(u>) approaches zero (for w k W[/2 and in some cases for other uj). One can showQ that near such a 
level, r(n) ~ \n — n cr | _1 ~ \t — t C r\~ 1 ^' 1 i-e. a sharp relaxation jump takes place accompanied by the generation of 
a burst of phonons (time t cr corresponds to n = n cr ). After a first jump, several other jumps may occur. For small 
n, formula (8) coincides with the result given by Fermi's Golden Rule (i.e. the standard perturbation theory). In the 
classical limit h — > 0, the decay of the mode energy vanishes for phonons in the ground state. 

The renormalized phonon density function p(u>) is nonzero also for lu larger than the maximal phonon frequency 
lom (due to phonon-assisted transitions described by the function f(t)), while p(to) is nonzero only in the one-phonon 
region < u> < u>m- Therefore formula (8) accounts not only for two-phonon decay of the local mode but also for its 
three- and more-phonon decay (the corresponding formula for a classical local modeQ accounts for two-phonon decay 
only). The linear and quadratic interactions considered above are commonly present, enhancing or weakening each 
other, as follows from the value of the resolvent R(u>) in (8). As a result, three- and more-phonon processes may also 
cause relaxation jumps, e.g. if loi > 2wm and two-phonon processes are forbidden by the energy conservation law. 
If u>i 3> 2lum, the consideration of linear and quadratic terms with respect to phonon operators in the interaction 
hamiltonian is not sufficient; one needs to include higher order terms (see, in this connection, ref.tj, where these terms 
have been considered within the classical description of the local mode). 
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To illustrate the effect of three- and more-phonon processes according to the quantum formula (8), we performed 
a model calculation of the relaxatipa rates of a Xe^-molecule in solid Xe. Note that the relaxation of this system was 
recently studied experimentally irJll where the relaxation jump was found to take place near the level n — 22. The 
level spacing in a Xej-molecule is equal to ^o;m(68 — n)/24. Two-phonon relaxation is allowed for upper levels until 
11 = 20. Lower levels can decay with emission of no less than three phonons (their spacing exceeds 2TwJm)- Only the 
symmetrized displacement of two nearest atoms along the Xe^ -molecule contributes significantly to the anharmonic 
interactions of this molecule with the environment (in this approximation, all matrices in (8) are scalars); the main 
contribution to this displacement comes from the longitudinal phoncpai3. The corresponding density function is well 
approximated by the cap-type shape p(x) = 8x^1 — x 2 , x — lu/lom^- Assuming f(t) — —Tt + j3(G(t) — G(0)), we 
calculated the r(n)-function for different values of the parameter (3 ~ Vu m (determining the shifts in the equilibrium 
positions of the phonon coordinates with the n — > n — 1 transition of the local mode) and with the maximum at 
n = 22. The results for 7(71) = T(n)uji n /ujM are presented in Fig.l. One sees that not only the two-phonon-decay but 
also the three phonon decay is strongly enhanced in the vicinity of the critical level n — 22. An estimation^ gives for 
Xe^-molecule in Xe crystal /3 w 0.2. 

Above the case T = was considered. Finite temperature effects are straightforward to take into account: in Eq. 
(8), there would appear an additional factor 1 + 2n(w) under the integral and the factor [e lu}it — 1) in f(t) should be 
replaced by the temperature-dependent factor (n(w) + l)(e lWit — l) + n(uj)(e~* UJit — 1) (here n(uj) — l/(e huJ / kT — 1)). The 
formulae obtained are applicable generally for nonadiabatic transitions in the case of linear diagonal and quadratic 
nondiagonal interactions of the system with the phonon continuum. 

In conclusion, we have presented a nonperturbative theory of multiphonon anharmonic quantum transitions between 
levels of a local mode taking into account linear and quadratic terms with respect to phonon operators. We showed 
that for high levels, transitions accelerate in time which is in contradiction with the standard perturbation theory. 
This continues down to the level n cr , when the rate of the process becomes large - of the order of , and a sharp 
relaxation jump takes place, accompanied by the generation of a burst of phonons. Depending on the level spacing, 
phonons may be emitted in pairs or triplets, quartets and etc. After the first jump, several other jumps may occur. 
Only for low levels standard perturbation theory and Fermi's Golden Rule are applicable. This theory can be applied 
generally for the nonperturbative description of nonadiabatic transitions caused by linear diagonal and quadratic 
nondiagonal interactions of the system with the phonon continuum. 

The author is grateful to Dr. M. Selg for his assistance concerning the enhancement of the three-phonon relaxation. 
This research was supported by Estonian Science Foundation, Grant No. 2274 and by the DAAD Grant HSPS. 



1 Ovchinnikov A. A., Zh. Eksp. Teor. Fiz. 57 (1969) 263. [Sov. Phys. JETP 30 (1970) 147]; Ovchinnikov A.A.and Erihman 
N.S., Usp. Fiz. Nauk 138 (1982) 290. [Sov. Phys. Usp. 25 (1982) 738]. 

2 Kosevich A.M.,Kovalev A.S., Zh. Eksp. Teor. Fiz. 67 (1974) 1793 [Sov. Phys. JETP, 40 (1974) 907]. 

3 Dolgov A.A., Fiz. Tverd. Tela (Leningrad) 28 (1986) 164. [Sov. Phys. Solid State 28 (1986) 907. 

4 Sievers A.J., Takeno S., Phys. Rev. Lett. 61 (1988) 970. 

5 Flach S., Willis C.R., Physics Report 295 (1998) 181. 
6 Klemens P., Phys. Rev.122 (1961) 443. 

7 Hizhnyakov V., Phys.Rev. B53 (1996) 13981; Proceedings of the Xllth Symposium on the Jahn-Teller Effect, Tartu 1994 
[Proc. Estonian Acad. Sci. Phys. Math. 44 (1995) 364]; Hizhnyakov V., Nevedrov D., Phys. Rev. B 56 (1997) R2908. 

8 Hizhnyakov V., Z.Phys. B 104 (1997) 675. 

9 Klinger M.I., Sov. Phys. -Usp. 30 (1987) 699. 

10 Lax M., Chem. Phys. 20 (1952) 1752. 

11 Hizhnyakov V., Selg M., Kink M., Kink R., and Maksimov J., J.Low Temp. Phys. Ill (1998) 709. 

12 Anselm A.L, Firsov Yu.A., JETP 28 (1955) 151; 30 (1956) 719. 

13 Lurie N.A., Shirane G., Skalyo J., Phys. Rev. B 9 (1974) 5300. 



4 



Figure caption 

Fig. 1. Theoretical relaxation rate 7 (in u>m ~ 8.3sec _1 units) vs level number n of a Xe^-molecule in solid Xe for different 
values of the phonon displacement parameter /3. 
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